%clear;

global C;
C = 1;

global omega;
omega = 1;

global eps;
eps = 0.1;

global u;
u = 1;

global n;
n = 2;

global q;
q = 1;

global h_;
h_ = 1;

global A;
A = inline('r', 'r');
%A = inline('r-r', 'r');

global inf_s;
inf_s = 2;

global h_s;
h_s = 0.5;

global inf_t;
inf_t = 2;

global h_t;
h_t = 0.5;

global inf_r;
inf_r = 2;

global h_r;
h_r = 0.5;

global h_theta;
h_theta = 0.5;

global h_alpha;
h_alpha = 0.5;

global inf_jt;
inf_jt = 1;

global h_jt;
h_jt = 0.5;

global phi0_cache;
phi0_cache = NaN(length(0 : h_theta : pi) + 1, length(0 : h_r : inf_r), length(0 : h_r : inf_r));

global isPhi0Cached;
isPhi0Cached = true;

global debugMode;
debugMode = false;
